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1 Introduction 



The Algebraic Bethe Ansatz (ABA) is a constructive method for solving 2d integrable 
models and it is natural to try to generalize the method to 3d integrable systems. The 
key relations in ABA are the Yang-Baxter equations (YBE), which ensure the commu- 
tativity of two transfer matrices with different spectral parameters. Zamolodchikov pre- 
sented a three-dimensional generalization of the YBE in IJ. His so-called tetrahedron 
equations (ZTE) imply the factorization of the scattering matrices of more than three 
one-dimensional objects - strings, thereby playing the same role for strings as the YBE do 
for the scattering of more than two particles in two-dimensional integrable field theories 
[21 E]- ZTE also generalize another function of YBE - integrability condition for the lattice 
models. They ensure the commutativity of two transfer matrices with different spectral 
parameters Qll^- In the three-dimensional case the scattering matrices (R-matrices) in 
the ZTE play the role of local evolution operators for the states on a 2d lattice (quantum 
spin model), or of the Boltzmann weights in the classical statistical models on 3d lattice. 
There exists a number of reformulations of ZTE in the context of integrable 3d models. In 
|31 the 3d conditions for integrability were formulated with the spin variables assigned 
to the lattice sites (interaction around the cube), while the spin variables in the original 
Zamolodchikov model were assigned to the faces of the 3d lattice. ZTE-equations for 
spin models where the spin variables live on the edges of a 3d cubic lattice (the so-called 
vertex formulation of ZTE) can be found in [71 EI- In the above listed articles a 
number of special solutions have been found and various parameterizations of ZTE were 
investigated (see also [13 [IHl [HI)- Finally 3d integrable models have been formulated in 
terms of free fermions in [H]. 

In the first section of this article we once more address the commutativity of two 
transfer matrices defined on a 3d lattice in the vertex formulation and we present a 
simplification of ZTE (STE), sufficient for the commutation of transfer matrices with 
different spectral parameters and defined on a 2d plane . The proof that STE implies the 
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commutation of the transfer matrices is similar the proof in the case of the full ZTE. We 
show that in a particular case the ZTE is formally equivalent to STE. In the last section 
we present a description of some solutions of STE which can be represented in terms of 
free fermions and we find the spectrum of the corresponding Hamiltonian. 

2 Commutativity of the transfer matrices. 

We consider a 3d classical statistical lattice vertex model on a cubic lattice with periodic 
boundary conditions. If -R-'^'f is the Boltzmann weight of a vertex which joins the links 
with variables labeled by i,j, k, and i,j, k, then the partition function Z of the model is 
defined as the product of the R's over all lattice sites: 

^ = E n ^ii (1) 

{i,j,k} all vertices 

We can rewrite it by defining a monodromy matrix T and the transfer matrix r, so 
that : 

^ = n'^y;t r = trr, Z = trr^^ (2) 

The monodromy matrix T acts on the N x N x N"^ dimensional spaces (we assume that 
the lattice is cubic with the equal numbers of sites in each direction: N^. = Ny = N~, = N) 
and Y\ denotes the product only over the horizontal (x y) plane of the 3d lattice. In the 
transfer matrix r the trace is taken over the N x N dimensional spaces of variables living 
on the edge sites of the horizontal (xy) lattice. The variables on the vertical links (labeled 
by Zij) can be regarded as an elements of the A^^ dimensional quantum space, while the 
horizontal ones form a N x N dimensional so-called auxiliary space, labeled by 
Fig. ^ for an illustration of the notation R^.y.z^^). 

In order to find a condition which ensures that two transfer matrices t{u) and t{v) with 
different spectral parameters, u and f , commute, we are going to define an extension of the 
YBE known from 2d models to the 3d models. In analogy with the method used for the 
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Figure 1: Rxi,yj,zij on 3D cubic lattice. 

YBE we will define two intertwiner R{u, v) matrices which will ensure the commutativity 
of the transfer matrices. At this point it should be mentioned that when we talk about 
spectral parameters in higher than two dimensions they should not be thought of as single 
real variable like in 2d, but as a set of continuous variables or angles on which the matrices 
Rijk depend. 

The transfer matrix constructed from the -Rij^-matrices can be represented as 

N N 
{xy} i=l i=l 

Let us denote the variables of the auxiliary spaces by and yj in the first transfer matrix 
t{u) and by a^, bj in the second transfer matrix t{v). As above the N x N states on the 
vertical links are denoted by Fig. Then the product of two transfer matrices is 

t{u)t{v) = (4) 

~ ^ ] iRxi,yi,ziiiu) ■ ■ ' Rxi,yff,Ziffiu)Rx2,yi,Z2ii'^)Rx2,yN,^2Ni^) ' ' ' "^a^iV.J/JV.^jViv ('^)) 
{x,y;a,b} 

X {Rai,bi,zii (v) ■ ■ ■ Rai,bN,ZNi ('^)-^a2,?ii,22i ('^) ' ' ' -Ra2,6]V,22JV ('^) ' ' ' RaM,bN,ZNN ('^)) ' 

We expect, that there is an intertwiner operator R{u,v) ensuring the commutativity of 
the transfer matrices with different spectral parameters: 



R{u,v)T{u)T{v) = T{v)T{u)R{u,v), 
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t{u)t{v) = tr {T{u)T{v)) = tr {R{u, v)T{v)T{u)R-\u, v)). (6) 

R matrix acts on the space labeled by the variables ai,bi,Xi,yi. Let us define the inter- 
twiner R matrix as the product of the local intertwiners: 

{u,v)Ry^^bj_{u,v) ■ ■ -Ryj^^bA^^v). (7) 

By use of a so-called railway argument one can prove the commutativity condition ®: 

tr {nu)T{v)} = tr v)T {u)T (v) R^^^J = (8) 

N N 

= Yl i^^N,aMiu,v)---Ry^,b^{u,v)Y[Y[{Rx,,y^,,^^{u)Ra,,y^,z,^{v)) 

{a,b,x,y} i=l j=l 

2 

{a,b,x,y} i=N 
N / N 

X Rx^,y^,zMRa^M,zM}WRy^,b^{'^^v) \ WRx,,y„z^,{u)Ra^^y^,z^^{v) 

1=2 \j=2 
N N \ 1 N 

1=2 j=l ) i=N i=l 

While the ordering of the Rxiaj and Rt.y^ operators is not important since they commute, 
we shall nevertheless treat them as non-commuting as long as possible (see later). 
If the following equations for local R-matrices are fulfilled for every i,j = 1, N: 

Rx,,aAu,v) Ry^,b,{u,v)Rx^,y^,z,,{u)Ra^^b„z,,{v) = 

= Ra,,b,,z,, {v)Rx,,y„z,, {u)Ry„b, («, v)Rx,,a, (m, v) . 

then by taking them into account in the figure brackets of the equation (jH)) and successively 
continuing that procedure along of all plane we will get 

2 

t{u)t{v) = Y YlRx,,a^{u,v){Ra,,buzAv)Rx,,yuzAu)Ry,M{u,v)Rx,,a,{u,v)} 
a,b;x,y i=N 

/ N 

X Ry^^b2{u,v)Rx,,y2,zi2{u)Rai,b2,zi2iv) iYlRy,,bA'^^'>j)Rxi,y.,z^Au)Ra^A,zu{v) 

\i=3 



X 



N N \ 1 Af 

i=2 j=l / i=N i=l 

2 / N \ ^ 

i=N \j=l J i=l 

N N \ 1 N 

^i=2 j=l / i=Af i=l 

= Yl^'^uh,,z,,{v)Rx^,y„z,,{u) J JJ i?:,.,,^, (^i, t^) (^a2,fei,^2l (^)^X2,j/i,^2i («)) 

\i=l / i=N 

N / N 

X Ryi,bA'^i'^)Rx2,y2(u,v) ,yj,Z2j \^)Ra2,bj,Z2j ip) 

i=2 \j=2 
N N \ IN 

X n n iu)Ra„y^,z,^ (v) R^.^a, (u, v) JJ R'^^^^ {u, v) JJ R'^^^ {u, v) 

i=3 j=l J i=N i=l 

(N N \ ^ ^ 

WWRa,,b,,z,,{.v)Rx^^y^^,^^{u) J WRy^^bXu,v) W Rx^^a,{u,v) 
i=l j=l J 1=1 i=N 

1 TV 

i=N i=l 

As Rxi,ai and -Ryi,;,^ commute it is seen from the last equation in (jlOp that 

t{u)t{v) = t{v)t{u). (11) 

Equations Q can be represented graphically as in We call them Semi-tetrahedron 
equations (STE) (referring to the known equations of Zamolodchikov). 

If we force the intertwiner matrices to have a structure similar to the Rijk (the Boltz- 
mann weights), we must consider the R'^ matrix elements as operators acting on an ad- 
ditional (auxiliary) space (the same space for all local intertwiners) : Rc,Xi,aii Rc,yi,bii and 
then take the trace in the transfer matrices over that space. In this case the intertwiner 
in the equation ((Tj) will be 

Ra,b;x,y;c 

{U,V) = Rc,xf^,yi^{u.,v) ■ ■ ■ Rc^xi,yi{u,v)Rc,aiMi'^^'^) ' ' ' Rc,af^ ,bfi{u,v) . (12) 
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Figure 2: Graphical representation of STE 

and the commutativity conditions of transfer matrices © are defined by Zamolodchikov's 
tetrahedron equations: 

Ra,,h,z,, {v)Rx,,y„z,, {u)Rc,y,A '^)Rc,x„a, {% v). 

Up to the last step, all equations in (jlUj) are valid in this general case as well, but the 
commutativity of intertwiners Rx„ai and Ry^^b, is broken for general Rc,xi,a, and Rc,yi,b, 
and we have to make a more detailed analysis. First we should take trace in the space of 
new auxiliary variable c 

tTc{Ry,b,c ® Rx,a,c}{Tx,y ® Ta,b} = {Ta,b ® T^,y} trc{Rx,a,c ® Ry,b,c}^ (l^) 

Rx,a,c Rc,Xi,yi 5 Ry,b,c Rc,Xi,ai • (-^^) 

i i 

As the variables x, a and y, b are independent, it is possible to change the positions of the 
R-matrices under the trace 

trc{-Rx,a,c ® Ry,b,c} = tTc{Ry,b,c ® Rx,a,c} (16) 



and following the by now familiar steps we arrive at the desired result 

tTc{Ry,b,c ® Rx,a,c}{Tx,y ® Ta^b}{i''^c{Rx,a,c ® -Ry,fe,c})"^ = ^^^^ 
= tTc{Ry,b,c <S> Rx,a,c}{Tx,y ® Tafi}{tTc{Ry,b,c ® Rx,a,c})~^ i 

tra,b,x,y{'^x,yiu) Ta,biv)} = 
= tT^a,b,x,y{tT^c{Ry,b,c ® Rx,a,c} {Tx,y{u) (g) Ta,b{v)}(tTc{Ry,b,c ® Rx,a,c})~^) = (^^) 

= tT^a,b,x,y{Ta,b{^) ® Tx,y{v)}. 

In the Figure 2 we present graphically the semi-tetrahedron and the tetrahedron equa- 
tions, respectively. It is seen that the semi-tetrahedron equations can be decomposed in 
two connected triangle equations. 

The semi-tetrahedron equations are cases of the tetrahedron equations. It is seen that 

if 

tv ci^RaxcRbyc) RxaRyb (-^^) 

the tetrahedron equations reduce to the semi-tetrahedron equations. Thus the structure 
of the STE is somewhat simpler than the structure of the ZTE and they consist of fewer 
equations than ZTE. 

3 Solutions of the STE using free fermionic represen- 
tations. 

We can present rather simple solutions of the STE by using a free fermionic realization 
of R-operator when the representation spaces where it acts are two-dimensional, i.e. the 
spaces of spin 1/2 states on the links of the 3d lattice. 

In the two-dimensional integrable models considered in ^21 ^1 ^| the three-particle 
Rijk matrices were represented in terms of fermionic operators: 

Ri23 : ® K2) O Ks) ® IJ2) ® Ijs), 
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Figure 3: Graphical representation of the semi-tetrahedron equations (the upper pair) 
and the tetrahedron equation (the lower pair). 



R 



123 



(^'^)(l,2,3) ~ ^' 



(20) 



with the fermionic creation and annihilation operators for each representation space de- 
fined as 

(|0)(0|), = l-n„ (|l)(l|), = n„ (|0)(1|), = Q, (|l)(0|), = c+, = c+c,. (21) 

Here p{a) = a are the parities of the states, due to the graded character of fermionic states 
[T^ . The R-matrices we will use here differ from the above matrices by a permutation 
operator 

Piall) ® |2) ® |3) = |3) ® |2) ® |1), (22) 
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i?i23 = P13R123 = i-iy^'^^^'^''^^'^''^^^^^^^^ (23) 

Similarly a representation for the two-particle R-matrix can be defined as 

R,, : |1)|2) ^ |2)|1), R,, = (-ir(^^)^(^-^ni?i2)l,t(bi)(^il)(b2)(^2|). (24) 

A representation of the three-particle R operators with /2jjfc-matrices which cannot be 
reduced to products of two-particle R-matrices was found in |i2j: 

i?i23 =: e^'"^-''^^'^''^ :, i,j = 1, 2, 3, an = 0, (25) 

where : () : means normal ordering. 

In general, when an 7^ 0, the representation above has 9 independent parameters aij 
(not counting the normalization parameter). On the other hand the -R123 operator in the 
so-called " free-fermionic" representation, : e^i,A°-^i~^ii^'^t'^3 (called so since the fermions 
enter quadratically) conserves the number of the fermions and has 20 matrix elements. 
So in this case there are 10 constrains between the matrix elements in fl23p . 

Similarly, the R12 operator : e^^^*-"'^"''*^'''^^'^^ : has 4 independent parameters a^j and 
one normalization parameter, so one constraint, known from the XX model is imposed 
on the six matrix elements 



pOOpll, plOpOl _ pOl plO (r)fi\ 

Let us write down all ten relations between the -R123 matrix elements, having the above 
" free-fermionic" representation. 

pOOl plOO I plOl pOOO _ plOO pOOl pOlO plOO I pllO pOOO _ plOO pOlO 

-'I'Oio-^'ioo "T -"'iio-'i'ooo ~ -'I'Oio-'i'ioo) -'T'ooi-^'ioo "T ^im^om ~ -'i-ooi-^-ioo) 

pOll pOOO I pOlO pOOl _ pOOl pOlO pllO pOOO I pOlO plOO _ plOO pOlO 

-'I'llO-'^OOO "T -"-lOO-^-OlO ~ -"-lOO-^-OlO' -"-Oll-^-OOO "T ^QOl^OlO — -"-OOl-^-OlO' 

pOll pOOO I pOlO pOOl _ pOOl pOlO plOO pOOl I plOl pOOO _ pOOl plOO 

-'I'lOl-^'OOO "T -"-lOO-^-OOl ~ -'I'lOO-'^OOl' -"-OlO-^-OOl "T -'^Oll-'^OOO ~ -'I'OlO-'^OOl' 

pllO pOOO I plOO pOlO _ pOlO plOO pOll pOOO I pOOl pOlO _ pOOl pOlO 

-'I'llO-^'OOO "T -"'lOO-'i'OlO ~ -'I'lOO-'i'OlO) -'T'Oll-^'OOO "T ^mi^QlO ~ -'I'OlO-^'OOl' 

plOl pOOO I pOOl plOO _ plOO pOOl 
-'I'lOl-'i'OOO "T -'T'OOl-'^lOO ~ -'T'OOl-'i'lOO' 
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pUUU Dili _ Jpj. 

-'^ooo-'^iii ~ 



DlOO 
-"-100 


dOIO 
-"-100 


p 


plOO 
-'I'OlO 


pOlO 
-'T'OlO 


p 


plOO 

-"■001 


pOlO 
-'T'OOl 


p 
K 



001 \ 

100 

001 
010 

001 , 

001 / 



(27) 



In terms of aij parameters equations contain 25 constraints - the independent ones 
from the 2^ x 2^ equations for matrix elements. Substituting ()23|) and (j^^ in the STE 
equations © we find 



0,1. (28) 



Note that when deahng with fermionic spaces as here, in the matrix representation of 
the vertex ZTE (I13p . some additional signs will appear in the equations (see the Appendix 
for details), compared to the usual form of these equations jH]. These signs are reflecting 
fermionic (graded) nature of our redefinitions of ZTE. 

Looking for solution for STE let us to choose the parameterization for R matrices as 



Rxa{u,V 



Rxyziu) =: e^j^-^^^^'^'^i :, 



. ^{'^]j{u,v)-5ij)c+ Cj 



(29) 



Then the fl28|) constraints can be rewritten as (we are omitting three equations, which are 
identities, and we are surpressing the arguments of the a*^ functions) 



'^21*^21 ~^ '^ll'^23'^31 



'^21*^21' 



1 3 2 3 4 

^21^12 ^ CL-y-^CL-y^CL 



32 ~ '^12*^21' 



2 3,2 3 4 
•^21*^22 ~r '^ll'^23'^32 



•^22*^21? 



13 13 4 

'^21*^11 + '^ll'^13%l 



1 4 
0210-11, 



nS'^SS ~ "'13"'22 "T "13"'33"'2l5 
2 4 

Qj-l O til 



"22"'12 ' "^^SU-S^"! 



''12"'12 



''23"'32"'ll5 



^23*^33 



3 2,4 3 2 
'^23*^22 + '^23%3'^21) 



1 4 

a,r,a 



3 2,4 3 2 
12^21 ~ '^2l'^12 + '^23%l'^ll5 



1 ^3 ^4 



'^22*^31 ~^ 



12"'33"'31 



'^33*^31' 



1 3 2 4 3 

'^22%2 + '^12%2%3 



'^33*^32' 



2 4 3 1,4 3 1 

'^12*^12 ~ '^12'^12 ' '^13'^32'^ll5 



2 4_3 2,43 2 
^12*^22 ~ '^22*^12 ~r '^23'^32'^11 ' 
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''12"'23"'32 



22^22 



*23"'32"'21' 



1 3 I 13 4 
'^22'^11 + '^12'^13'^31 

12 13 4 

'^22*^21 "I" '^12'^23'^31 



2 3,2 3 4 
'^22*^12 ~r '^12^13*^32 



21^32 



2 3 4 
+ 011033032 



^32) 



3 1 4 3 1 

'^ll'^22 + '^13'^3l'^2l5 



3 2,4 3 2 

— '^2l'^22 ~r '^23'^3l'^21' 



^12^22 "I" ^13^32^21) 



13 13 4 

O21O31 + OiiOs^aci 



''11"'33"'31 



^31' 



'^23 — '^23*^12 "I" '^9sO 



23"'33'^11' 



'^13 — '^13*^12 "I" '^13'^ 



13"'33'^11' 



4 2 
a2iO,ii, 



2 4 

•^11^12 



4 1 
012011- 



In the general case one has the following constraints on the a^j parameters 



ni 



(3) (3) 
^23 '^32 _ 

"13 "31 

,(3) 



*23 



(2) 

— Oil , 


^22 = a 


(2) 
22 ) 




(4) (4) 
'23 '^32 


(3) (3) 
^23 '^32 


(3) 
'^22 


(4) 
•^22 


(4) (4)' 
'13 %1 


"13 "31 


«f,' - 


"11 










"23 




«!? ' 





"11 — 



"21 — 



ni 



ni 



a(4)a(3) _yl3(4)' 
"33 "11 A13 



,(2) 



(3) (4) 
%1 '^13 



J4) (3) _ 13(4)' 
"33 "11 A13 

(3) (4) (4) 
(2) ^ Q23 Q3I ^13 

" ag'(4'<.f,'-xl3'")' 

^fcr(«) _ „(")„(«) 



^22 ~ 



"12 — 



(2) 
«21 = 



„{4) 13(3) 
"33 A 13 



„(3) 13(4) 
"33 A 13 



a(4)a{3) _ yl3(4) 
"33 "11 A 13 



(4) (3) 
'^31 '^13 



O33 



^13(4) ' 
A13 



(4) (3) (3) 
'^23 %1 '^13 



^23^ ('^33'' '^11'' X13*' "*) 



(«) (a) 



3,4. 



We will pay special attention to the homogeneous case 



(r) _ (r) 



r = 1, 2, 3, 4, = a>2^, r = 1, 2, a^i^ = 0^2^ = , a = 3, 4, 

where the above equations can be converting into the following equations: 

= const = A, q; = 1, 2, r = 3, 4, 



12(0!) , 1 

xli + 1 



(r) 



,(") _ „(«) 



(a) 



a.: 



a. 



(r) 
33 

(2) ^(r) 



^3 



13 



(a) 
"11 

,(0 
''23 ' 



(r) 
«12 



aji^ + A, r = 3, 4, A = const. 
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The solution of these equations for R2 can be parameterized via trigonometric functions 
Oil = 0^22 = sm{u)/ cos(u), ai2 = 021 = 1/ cos(u), A = 0. (35) 
This is the parameterization of the XX-model. For for the R3 we then obtain 

ai2 = 0,21 = an = 0-22 = ^33 = sm{u)/ cos(m), 

ai3 = a-ii = = ^32 = 1/ cos(m), a = 0, (36) 

which is an extension of the XX R-matrix to three dimensional case. Then, by using the 
technique of coherent states developed in jHl El one can find a fermionic expressions of the 
transfer matrix t{u) and the partition function Z = n^''"('^)- 



TiU 



JV JV 

/2 ' 2 
Dipntfj : exp{^ [ai2i>x,i,ji^y,i-i,j-i + a2i^y,i,ji)x,i-i,j-i + aiiV'a;,ij^x,i-2,: 



i,jCz,i,j-l 

+ a32Cz,i,j'^y,i,j-l + ('^33 " ^)Cz,i,jCz,i,j - ^x,i,j^x,i,j " i^y,i,j'^y,i,j)} '■ ■ (37) 

After Fourier transformation we will have 

JV JV 

/2 ' 2 
DipDip Y[ ■ exp{(ai2e*(^+^^V'^,p,g?/^y,p,g + a2ie*^^+''^V'j;,p,g^x,p,q + aiie^'^^px,p,q4^x,p 

~l~ 0,22^ '^'4^y,p,q'^y,p,q ~l~ ai3€^1px,p,qCz,p,q ~l~ O^-iC^ Cz^p^q'ipx,p,q Oj23€ '^Ipy^p^qCz^p^q 

~\~ Oi32(^ Cz^p^qlpy^p^q ~\~ (0-33 ^")(^z,p,qCz,p,q '4^x,p,q'4^x,p,q '^y,P,q'^y,P,q^} 

JV JV 
2 ' 2 

= n "'^P' ^) ■ exp{v4p,gC2,p,5C2,p,5} : . (3 
Here the normal ordering notion :: refers to the Cz and Cz operators and 

a{p, q) = 1 + (aiia22 - 012021)6^*^^+''^ - (ane^'^ + 0226^'^), 

X 

Ap,q = 033 - 1 - ((0226^''^ - l)ai3a3ie2*P + {aue^'P - 1)0320236^*" 

- ai2a3ia23e2^(^+^) - a2iai3a32e''^P+'^^) . (39) 
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The solution corresponding to (j36p reads 

a{p,q) = l-sm{u)/cos{u){e^'P + 6^''^), (40) 

A„a = sm(u)/ cos(m) — 1 H r^rr, r-r-, , w n- ^-tt- 

^'^ ^ " ^ ' cos(m)2(1 -sin(M)/cos(M)(e2*P + e2*9)) 

At the point m = we have 

a{p,q) = l-u{e^'P + e^"^), 

Ap,, = u-l + {e^'P + e^"^){l + u{e^'P + e^"^)), (41) 

and the transfer matrix r{u) takes the following form for small values of the spectral 
parameter u 

JV JV 
2 ' 2 

JV JV 
2 ' 2 

JV JV 
2 ' 2 

— g-P(P,g)5z,p,9C2,p,5+M£(p,<j)c2,p,qC2,p,q ^^2^ 



where 



p,q 



^nP,i) = gi^P + gi^g^ ^(p^ = 2 cosh P{p, q). (43) 

When — — oo ( or — — oo) this transfer matrix decouples into the product of N 
independent XX model's transfer matrices defined on a ID chain. 

In the case an = there are another solutions for the fermionic STE equations (j3T| . 
They are defined by three multiplicative relations for the parameters of three-particle 
R-matrices 

4^a'l2 = (3, a^3 = aa^2a^3, flgi = 7a32«2i, k = 3,4, (44) 
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where the constants a, (3, 7 can be considered as model parameters. The relations can be 
reformulated as 

a2i{u)ai2{u) = /3, ai3{u) = aai2{u)a23{u), a3i{u) = 7032 H 021 H- (45) 

The non zero matrix elements of corresponding Rijk operators in this case are the 
following (with above relations) 



dOOO 
-'T'OOO 


= 1, 


pOOl 
-'^100 


= 012, 


plOO 
-n-OOl 


= O21, 


pOll 
-'T'Oll 


= 032023, 


pOlO 
-'T'lOO 


= Oi3, 


plOO 
-n-OlO 


= O31, 


pllO 
^110 


= 013031, 


pllO 

-"-oil 


= — 023031, 


pOll 
-n-llO 


— "013032 


plOl 
-'^101 


= —012021, 


pOll 
-"-101 


= 023012, 


plOl 
-"-110 


— ~03iOi2 


plOl 
-n-Oll 


= —032021, 


pOOl 
-'T'OlO 


= 032, 


pOlO 
-n-OOl 


= -O23, 


pllO 
-n-lOl 


= —013021, 


pill 
-n-iii 


— 013032021 + 023031012. 







(46) 



While the non vanishing matrix elements of two {R^''^)ij intertwiners as solutions of the 
semi-tetrahedron equations have the following expressions in terms of the Uij functions of 
R123S 

{R^{u,v))'^^ = 032(M)/o32(t;), {R^{u,v))\\ = a2i{v)/a2i{u), 

{R\u,v))ll = 032(^)021 (^;)/o32(^)o2i(k), {R\u,v))% = 1, 

(47) 

(i?2(ti,^;))00 = a^2(u)/a,2(v), (R''{u,v))l\ = 023(^;)/023(^i), 

{R'^{u,v))l1 = Ol2(^^)o23('y)/ol2(^;)o23(^^), {R^{u,v))% = 1. 

As in two dimensional case we can extract from the commuting transfer matrices 
exact solvable quantum models (in our case on a two-dimensional lattice rather than on 
one-dimensional lattice) , identifying first the logarithmic derivative of the transfer matrix 
with the Hamiltonian: 

\nr{u) = In i tr^^^y. JJ R^^^y^^,^ .{u) J = 

V id / 

= lnr{uo) + u{r{uo))~'^dr{u)/du\u=uo + ■■■■ 
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By inserting the representation for R^.^y-^^i in this equation and taking the 
trace over the auxihary spaces, due to the gaussian form of the exponents in R-operators, 
r will have the similar form (in the following we shall omit the Zij notation, keeping only 

ij) 

ij r i,j r i,j r 

= {PY~\- + -)ana3i, Br = {PYaisan, = (/9)^a32a23- (48) 

a 7 

From (22) it follows that 

023032 = 013031/(0/37). 

If 013031 is small it is possible to expand r in terms of the 013031 variables and H in 
the expression above manifests itself as a non-local Hamiltonian. From the structure 
of H it follows that although Hamiltonian is not local, there is localization around the 
path (z, j — s> i + 1, j + 1 — s> ... — s> i + r,j — r...). In momentum space this becomes more 
apparent 

'^^^P^,Py' ^P^'Py) = ^''a ~^ 7 ~^ ~^ 07^ ^'')'^i,P:.+Po' '^P=o,P.+PO^ 

PT,Py=Px+PO 

Po = t\nf3. (49) 

In the case where 

rv 

a7 = (/3)^ - = (50) 

7 

the Hamiltonian is real and can be written as 

'^(4.,Py^ ^V..Py) = 2 ^(Coshpi + COSp)c+p+p„, Cp^p+p^, 

P 

01 

Po = In 07/2, pi = ln— . (51) 

7 

In the limits /? — > 0, 00 Hamiltonian PHj) becomes completely local. 
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The expression (15 1|) allows us to represent the transfer matrix as r = expTi, and 
thus to find the exact discrete time Hamiltonian not only for small values of 012021- The 
Hamiltonian follows directly from: 

T = A:e":=AY[: e^W^^'= = Ae^W^-^fe, e{k) = In (1 + e{k)). 

k 

As one can see from the expression ()51|) rotational symmetry is broken in the above 
model. In order to recover it we can suggest another integrable model , where the transfer 
matrix is constructed with the help of the same solutions -R3 and R2 (satisfying (|46p and 
dUj)) of the local STE equations Q, as in the previous section, but with some chess like 
disposition of the R3 and rotated -R3 matrices 

f{u) = trf{u) = 5^ n ( n^-«-^^.--.^.(")^"..+....^..+i,j^) 

{x,y} j \ i 

n -^a;2i,y2j+l,Z2i,2j+l(^)-^a;2i+i,i/2j+l,2;2j+i,2i+i('")^ ' (^^) 

where r^, i = 1, 2, 3, refer to the rotated matrices 

TDTii'j'k' _ Tji'jk' 
ijk ~ ^ij'k ' 

^"S'' = (53) 

UT3i'j'k' _ Tjijk 
ijk ~ ^i'j'k- 

In this case the multiplicative intertwiner matrix, necessary for the commutativity of two 
transfer matrices with different spectral parameters: 

Ra,b;x,y{u, v)f{u)f{v) = f {v)f (u) Ra,b;x,y{u, v) , (54) 

is modified a little compared to ((Tj) 

j i 

Here r denotes 

= Rfa (56) 
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The commutativity can be verified by putting the right hand sides of (j52p and (j55|l 
into the equation and observing that it is fulfilled if four local equations are satisfied 
(as it was done in (jHJ and (fTUI) ). Here we are omitting all arguments of the operators: 

RybRxaRxyzRabz RabzRxyzRxaRyb^ l^"^) 

/ P TxTa\—lT) p ''"^^ P '''a P '''a P '''^ f? ( TiTaN— 1 ( 

y-l^xa J J^yb^xyz ^abz J^abz J^xyz J^yby^^xa ) i v"-'"/ 

i^Ryb ^ ) RxaRxyz ^ Rabz Rabz Rxyz ^ Rxai^Ryb ^ ) > (^^) 

{Ryb^^^T\RxJ''^T'Rxyz^^^'Rabz^''^' = Rabz^''^'Rxyz^^^'{Rxa^''^T\Ryb^'^T'- (60) 



By Tx,Ty,Ta, Tb we mean the rotations or transpositions defined in and 



p TxTa _ 

^xa 


_ P T 

^xa 1 


Ryb ^ Ryb 5 




p Tx 

^xyz 


R ^2 

^ ^xyz ) 


P Ta p T2 

'^abz ^abz 1 




^ ^xyz 


T3 Tl 

J^xyz 1 


P Tfc p Tl 


(61) 


P TxTy _ 

^xyz 


^xyz ) 


T3 TaTb p T3 

-'■ '^abz -'■ '^abz 





The solutions of the first equation (fS7j) (which are the STE defined in (P)) are also 
solutions for other three equations. It can be shown for ()58p by applying corresponding 
r operation on the standard STE equations as done below: 

(^RybRxaRxyzRabz) (RabzRxyzRxaRyb) 1 (62) 

T? R R R '''xTa p TxTa p '''a f? f? ( f{'^\ 

^yb-'^xyz ^abz ^xa ^xa ^abz ^xyz -'^ybi \^'~' ) 

Two other equations, and (jUUI) . can be derived analogously. 

The corresponding 2d quantum Hamiltonian in the fermionic representation for this 
model, unlike the previous one, has only nearest-neighbor and next to nearest-neighbor 
interactions, and exhibits a Manhattan like structure on the 2d square lattice shown in 
Fig. El 

f = (l-/?2)^:eT^^('="^^, (64) 
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+ a23a32 ^(c^i,2j+lC2i,2j + /^Q;C^i-l,2i+lC2j,2j + /^C^i-l,2jC2i,2i + P'^ Oi<4i,2j'^2i,2j) + • 

The H{c'^,c) consists of mass terms and hopping terms only along the elementary 
circuits, described by the arrows drawn on the links of the square lattice in the Fig|21 
(if continued for all the links, they form a square Manhattan lattice). In we find 
terms with C2i,2j operators. They correspond to the circulations around two elementary 
circuits with clockwise and counterclockwise orientations, connected by the 2i, 2j vertex. 
The other three terms with the annihilation operators on the even-odd, odd-odd, even- 
even sites can be obtained from that expression, following to the directions of the arrows, 
outgoing from that site and writing the same horizontal, diagonal and vertical hopping 
parameters (and also the mass terms) as in the first term. 

In the momentum space 

2jp+2kq _l_ 1 AT/O r,_2jp+2kq I 



(^2j,2k — l^p,q=l ^ ^lp>g' ^2j,2k — N ^P,q=l ^ ^^P,q^ 

C2j+l,2k — Z^p,q=l 6 ^ '^2p,q, ^2j+l,2k ~ N l^p,q=l ^ 



Co, 



C2j+l,2fc+l - ]v 2^p,q=l ^ '^3p,g, C^^+l,2k+\ " TV 2^P,g=l ^ 

,+ 



_ 1 ^iV/2 27r52£±i^±ilS„ „+ _ 1 sr^l2 2n22I+^h+lll + 

C2j,2k+1 - jv 2^p,q=l ^ ^ '^4p,5, %-,2A:+l " JV l^p,q=\ ^ ^ C, 



(65) 



the Hamiltonian ()65|) will take a more compact form 

E{j), q) = c^p gC2p,g(a3iai3/5e"'P + 0230326*^) + 4p,gCip,g(a3iai3e*^ + 023^32/56"*^) 

+ 4p,qC3p,q{a3iai3^e''P^''' + a23a32/3ae'^-*') + 4^^^01^,5(031013^6'^"*'' + a23a32/5ae"*^+''') 



xq\ 



+ ctpqC4p^g{a-iiai3pe"^ + 033032/36 *") + 4^ gCip,g(a3iai3/3e + 0230326 

+ 4p,gC3p,g(03iai3/5e*^ + 0230326"*'^) + 4p,gC2p,g(a3iai3e~*^ + 023032/36*^) 

+ 4p qCAp,q{a3iai3(3ae'P+"^ + 023032 -e"*^'"'^) + ct „C2p,q{a3iai3(3ae~'P~"' + 023^32 -e'^+*'') 

a a 

+ 4p,gC4p,g(a3iai3/5e*'' + 0230326"*'') + 4p,5C3p,g(a3iai3e"*'' + 023^32/36*'') 

+ (4p,gClp,g + 4p,<7C3p,g)(a3iai3^ + 023032/?^") (66) 
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Figure 4: Manhattan like structure on 2d square lattice 

+ {'^tp,q'^2p,q + clp,qC4p,q){o.3lO'13P'^Oi + 023032 — ). 

For the simple choice of the parameters 031013 = O23O32 — a^, a — 1, (3 
Hamiltonian H[p, q) acquires simple matrix form 

/ 1 cosp cos (p — q) cosq ^ 

cosp 1 cos q cos {p + q) 

cos {p — q) cos q 1 cosp 

y cosg cos(p + g) cosp 1 J 

two nonzero eigenvalues of which gives us the spectrum 

e± = 80^(1 ± cosp cos g). 



1, the 



2o^ 



(67) 



(68) 
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Appendix: ZTE equations for free fermionic representation ( 1291) 
of R- operators 

In this appendix we present the independent set of ZTE equations (fT^ in case when 
the three state R-operators admit the : e^°'^^~^^^''^^^^ : free-fermionic reahzation. For the 
operators (j23|) the equations (jl3j) . written in matrix form, are: 



In the equations we have omitted spectral parameters and we have distinguished the 
elements of different matrices by labels 1, 2, 3, 4. In terms of a^j elements the non-identity 
equations are 




/ V ^^a^b^-z ^^xi'yjz ^cjyjb ^Jcjxja ^ ^ 




(69) 




(70) 



'SS'^Sl + '^13'^2l'^32'^33 + '^13'^3l'^32 + '^23'^3l'^33 ~ %l'^33' 
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^13^11 ~ ^21^13^11 "I" ^31^23^12^11 + ^31^13^12 + ^11^13' 
•^22*^12 ~ '^12'^22 ~^ '^22'^12'^21; 
^22^32^13 ~^ ^32^12 ~ '^22 '^12 '^31 "I" ^12*^32) 
'^22*^32^23 "I" '^32*^22 ~ '^22'^32' 

^22^32^33 "I" ^32^32 ~ ^32) 

2 3 1 3 2 1 

^12 ~ '^i2'^12 + '5'22'^'12'^'ll! 

2 3 421 41 4321 431 

•^23*^12 ~ '^22'^13'^21 ^" '^12*^23 ^" '^32'^23'^12'^21 ^" '^32'^i3'^22' 

'^23^32'^13 "I" '^33'^12 ~ '^22'^13'^31 + '^12'^33 + '^32'^23'^12'^31 + ^32^13^32' 

'^23*^22 ~ '^22*^23 "I" ^32'^23'^22) 

^23^32^23 ~ ^32^23^32' 

^23*^32*^33 + '^33*^32 ~ '^32*^33' 
•^13 ~ '^22'^13'^11 + '^12'^13 + '^32'^23'^12'^11 + '^32'^13'^i2 ) 



(71) 



'^13 ~ '^23'^13'^21 + ^13'^23 + "^33 '^23 '^12 '^21 + '^33 '^13 '^22 ) 
^33^^13 ~ ^23^13^31 "I" ^13^33 "I" '^33'^23^12'^31 "I" ^^33^^13^32' 



^23 ~ ^3^23^22 + ^23^23' 

^33*^23 ~ '^33*^23*^32 + '^23*^33' 
~ '^23'^13'^11 + '^13'^13 ~^ ^33'^23'^12'^11 + '^33'^13'^12' 

^Sl'^ll "I" '^ll'^21^32'^13 "I" ^^11^^31^12 + '^2l'^3l'^13 ~ ^11^31' 
'^2l'^21 "I" '^ll'^2l'^22 ~ ^21) 



O.^l Ooi + all ttoi Cloaaoq + C ji Coi + Ool O'?! floQ 



''31"'21 "'11"'21"'32"'23 "'11"'31"'22 "'21"'31"'23 ~ '^31' 
41*^31 + '^ll'^2l'^32'^33 + '^ll'^3l'^32 + '^2l'^3l'^33 ~ 0- 
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